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We present the new analysis of the semileptonic B decays in the framework of the 
relativistic quark model based on the quasipotential approach. Decays both to heavy 
DW and light vr(p) mesons are considered. All relativistic effects are systematically 
taken into account including contributions of the negative-energy states and the wave 
function transformation from the rest to moving reference frame. For heavy-to-heavy 
transitions the heavy quark expansion is applied. Leading and subleading Isgur-Wise 
functions are determined as the overlap integrals of initial and final meson wave 
functions. For heavy-to- light transitions the explicit relativistic expressions are used 
to determine the dependence of the form factors on the momentum transfer squared. 
Such treatment significantly reduces theoretical uncertainties and increases reliability 
of obtained predictions. All results for form factors, partial and total decay rates 
agree well with recent experimental data and unquenched lattice calculations. From 
this comparison we find the following values of the Cabibbo-Kobayashi-Maskawa 
matrix elements: \V cb \ = (3.85±0.15±0.20) x 1(T 2 and \V ub \ = (3.82 ±0.20 ±0.20) x 
10~ 3 , where the first error is experimental and the second one is theoretical. 

PACS numbers: 12.39.Ki, 13.20.He, 12.15.Hh 



I. INTRODUCTION 



The investigation of the semileptonic decays of heavy B mesons is an important source 
for the determination of the parameters of the standard model, such as Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements V cb and V u f,- They also provide valuable insight in quark 
dynamics in the nonperturbative domain of QCD. Recently significant experimental progress 
has been achieved in studying exclusive semileptonic B decays (for a recent review see, e.g., 
Ref. [1]). A very important information on both the values of the total decay rates and 
the differential decay rate dependence on the momentum transfer is becoming available. 
The experimental accuracy is constantly increasing due to the large data accumulated at B 
factories. For most decay modes experimental errors are comparable to or already smaller 
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than theoretical ones. Thus the accurate extraction of the CKM matrix elements from the 
semileptonic B decays requires an increase in the reliability and precision of theoretical 
methods for determining weak decay form factors. Main theoretical approaches for calcu- 
lating these form factors are lattice QCD, QCD sum rules and constituent quark models. 
Unfortunately, the first two approaches are applicable only in limited regions: lattice QCD 
gives reliable results for high values of the square of momentum transfer from the parent to 
daughter hadron q 2 , while QCD sum rules are suitable only for low q 2 . Thus extrapolations 
for the full q 2 region are at present inevitable in these methods. On the other hand, quark 
models, which are not straightforwardly related to QCD (at least this relation is not cur- 
rently established), allow to determine weak form factors in the whole kinematical range. 
To give a correct description of semileptonic decays such models should consistently and 
comprehensively account for relativistic effects and use the wave functions which lead to 
correct meson masses. They should also respect all the relations imposed by the symmetries 
of the QCD Lagrangian arising in the heavy quark limit. 

From a theoretical point of view the simplest semileptonic B decays are heavy-to-heavy 
transitions such as B — > D^eis. The presence of the heavy quark both in the initial and 
final mesons significantly simplifies the understanding of such processes. The heavy quark 
limit rriQ — > oo is a good initial approximation [2]. In this limit heavy quark symmetry 
emerges substantially reducing the number of independent characteristics which are neces- 
sary for the description of heavy-light meson properties [3]. Mass and spin decouple from the 
consideration and all meson properties are determined by the light quark degrees of freedom 
alone. This leads to symmetry relations between form factors responsible for the heavy- 
to-heavy weak transitions. Thus one needs considerably less independent functions. For 
the semileptonic B decay to the ground state meson all form factors can be expressed 
through one Isgur-Wise function [2] which is normalized to 1 at the point of zero recoil of 
the final meson. However, in reality b and c quarks are not infinitely heavy and therefore 
the corrections in inverse powers of the heavy quark mass uiq (especially m c ) are important. 
Heavy quark effective theory (HQET) [4] is the adopted tool for a systematic expansion 
of weak decay amplitudes in I/tuq. The coefficients in this expansion are functions of the 
velocity transfer in the weak decay and do not depend on spin and flavour of the heavy 
quarks. These functions originate from the infrared (nonperturbative) region of QCD and 
thus cannot be determined from first principles of QCD at present. It is necessary to use 
some model assumptions in order to calculate these functions. Note that notwithstanding 
heavy quark symmetry violation already in the first order in 1/m.Q, its remnants remain 
and thus HQET significantly restricts the structure of I/tuq corrections. E.g., in the first 
order in 1/toq for B decays to the ground state mesons one mass parameter and four 
additional functions (two of which are normalized at the zero recoil point) emerge instead 
of 12 possible ones. This is the consequence of QCD and heavy quark symmetry. Thus all 
models which pretend to describe correctly weak heavy-to-heavy transitions should satisfy 
HQET symmetry relations. 

The theoretical description of exclusive semileptonic heavy-to-light B decays, such as 
B — > Tr(p)ev, represents a more difficult task since the final meson contains light quarks 
only. The expansion in inverse powers of the b quark mass does not reduce the number of 
independent form factors. Only relations between semileptonic and rare radiative decays 
emerge in the heavy b quark limit. It is important to note that in these decays the final 
light meson has a large (compared to its mass) recoil momentum (energy) in the rest frame 
of the decaying B meson almost in the whole kinematical range except the small region near 
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the point of zero recoil. The maximum value of recoil momentum is of order of mj/ '2. Thus 
near this point one can expand both in inverse powers of the heavy b quark mass and large 
recoil momentum of the final light meson. Such expansions lead to the so-called large-energy 
effective theory (LEET) [5] , to a new symmetry and as a result to the form factor relations 
in the heavy quark and large recoil limits [6, 7]. Large values of the recoil momentum require 
the completely relativistic treatment of the semileptonic heavy-to-light B decays. 

In our previous papers [8-10] we considered semileptonic B decays to heavy and light 
7r(p) mesons in the framework of the relativistic quark model based on the quasipotential 
approach. At that time taking into account large experimental errors in the measured decay 
rates we used the simple Gaussian parameterization for the meson wave functions. Moreover, 
we calculated the heavy-to-light decay form factors only at the point of maximum recoil of 
the final light meson using an expansion both in inverse powers of the heavy b quark and large 
recoil momentum of the light meson. Then we employed a Gaussian or pole parameterization 
of the form factors in order to extrapolate them to the whole kinematical range. Such 
substitutions and extrapolations induced significant theoretical errors in the obtained results, 
but the accuracy was sufficient compared to large experimental uncertainties. Since then, 
as it has been already mentioned before, the experimental accuracy improved significantly. 
The main aim of this paper is to revise our previous considerations of semileptonic B decays 
substantially increasing the precision and reliability. This is achieved by using the wave 
functions of the heavy B, D and light ir, p mesons which were obtained by calculating their 
mass spectra [11, 12]. The complete expressions for the heavy-to-light decay form factors 
are used for calculations and the determination of their q 2 dependence, thus avoiding ad hoc 
parameterizations. In the following we concentrate on the study of the relativistic effects 
and, for simplicity, neglect short- distance radiative corrections [4] since their contribution 
does not exceed the uncertainty of our calculations. 

The paper is organized as follows. In Sec. II we briefly describe our relativistic quark 
model. Then in Sec. Ill we discuss the relativistic calculation of the decay matrix element of 
the weak current between meson states in the quasipotential approach. Special attention is 
devoted to the contributions of the negative energy states and the relativistic transformation 
of the wave functions from the rest to the moving reference frame. Semileptonic B decays to 

£>(*) 

mesons are considered in Sec. IV using the heavy quark expansion in 1/rriQ. Leading 
and subleading Isgur-Wise functions are explicitly determined as the overlap integrals of the 
initial and final meson wave functions. A comprehensive comparison with recent experimen- 
tal data is given and on this basis the value of the CKM matrix element \V c b\ is determined. 
In Sec. V semileptonic B decays to rr and p are investigated. The parameterization of the 
calculated form factors in the whole kinematical range is given. Total and partial decay 
rates are compared with recent measurements and the value of the CKM matrix element 
\V u b\ is extracted. Section VI contains our conclusions. 



In the quasipotential approach a meson is described by the wave function of the bound 
quark- ant iquark state, which satisfies the quasipotential equation [13] of the Schrodinger 



II. RELATIVISTIC QUARK MODEL 



type [14] 
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where the relativistic reduced mass is 



EiE 2 M 4 - {ml - m 2 ) 2 

and Ei, E 2 are the center of mass energies on mass shell given by 

M 2 — m\ + m\ M 2 — m\ + m\ 

1= 2M ' 2= 2M ' ^ 

Here M = E\ + E 2 is the meson mass, m\ y2 are the quark masses, and p is their relative 
momentum. In the center of mass system the relative momentum squared on mass shell 
r Gelds 

h 2 (M] = W 2 ~ K + m 2 ) 2 ][M 2 - (mi - m 2 ) 2 ] 

4M 2 • ^ 

The kernel V(p, q; M) in Eq. (1) is the quasipotential operator of the quark- ant iquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by [15] 

V(p, q; M) = iii(p)«2(-p)V(p, q; M) Ul (q)u 2 (-q), (5) 

with 

V(p, q; M) = ^ s ZV(k) 7 ^ + Konf (Wr 2;M + V c ^(k), 
where a s is the QCD coupling constant, is the gluon propagator in the Coulomb gauge 

D oo (k) = _g ) D ii^) = -^(&-!™y D« = £>» = 0, (6) 
and k = p — q; 7 M and u{p) are the Dirac matrices and spinors 

u x (p) 



> e[P) + m \ °P I X A - (7) 

Here cr and x X are the Pauli matrices and spinors; e(p) = \/p 2 + m 2 . The effective long-range 
vector vertex is given by 

r M (k) = 7^ + ^<v*"» (8) 

where k is the Pauli interaction constant characterizing the long-range anomalous chromo- 
magnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 

V v (r) = (l-e)(Ar + B), 

V s (r) = e(Ar + B), (9) 
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reproducing 

V CODl (r) = V s (r) + V v (r)=Ar + B, (10) 

where e is the mixing coefficient. 

The expression for the quasipotential of the heavy quarkonia, expanded in v 2 /c 2 without 
and with retardation corrections to the confining potential, can be found in Refs. [15, 16], 
respectively. The structure of the spin-dependent interaction is in agreement with the param- 
eterization of Eichten and Feinberg [17]. The quasipotential for the heavy quark interaction 
with a light antiquark without employing the nonrelativistic (v/c) expansion for the light 
quark is given in Ref. [11]. All the parameters of our model like quark masses, parameters of 
the linear confining potential A and B, mixing coefficient e and anomalous chromomagnetic 
quark moment k are fixed from the analysis of heavy quarkonium masses [15] and radiative 
decays [18]. The quark masses = 4.88 GeV, m c = 1.55 GeV, m u ^ = 0.33 GeV and 
the parameters of the linear potential A = 0.18 GeV 2 and B = —0.30 GeV have usual 
values of quark models. The value of the mixing coefficient of vector and scalar confining 
potentials e = — 1 has been determined from the consideration of the heavy quark expansion 
for the semileptonic B — > D decays [8] and charmonium radiative decays [18]. Finally, the 
universal Pauli interaction constant k — —1 has been fixed from the analysis of the fine 
splitting of heavy quarkonia 3 Pj- states [15]. Note that the long-range magnetic contribu- 
tion to the potential in our model is proportional to (1 + k) and thus vanishes for the chosen 
value of k — — 1. It has been known for a long time that the correct reproduction of the 
spin-dependent part of the quark- antiquark interaction requires either assuming the scalar 
confinement or equivalently introducing the Pauli interaction with k — — 1 [15, 16, 19] in 
the vector confinement. 



III. MATRIX ELEMENTS OF THE WEAK CURRENT FOR b -» c, u 

TRANSITIONS 

In order to calculate the exclusive semileptonic decay rate of the B meson, it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
In the quasipotential approach, the matrix element of the weak current = /7 M (1 — 75) ft, 
associated with b — > / (/ = c or u) transition, between a B meson with mass Mb and 
momentum ps and a final meson F (F = or 7r(p)) with mass Mp and momentum pp 
takes the form [20] 

{F(p F )\jJ\B(p B )) = J ^^^ PF (p)r M (p,q)* Bps (q), (11) 

where T M (p,q) is the two-particle vertex function and $m Pm are the meson (M = B,F) 
wave functions projected onto the positive energy states of quarks and boosted to the moving 
reference frame with momentum p«. 

The contributions to T come from Figs. 1 and 2. The leading order vertex function 
corresponds to the impulse approximation, while the vertex function T^ 2 ^ accounts for 
contributions of the negative-energy states. Note that the form of the relativistic corrections 
resulting from the vertex function T^ 2 ) is explicitly dependent on the Lorentz structure of 
the quark- antiquark interaction. In the leading order of the the heavy quark expansion 
(^fe, c — ► 00) for B — > D transitions only contributes, while contributes already at 
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FIG. 1: Lowest order vertex function contributing to the current matrix element (11). 
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FIG. 2: Vertex function taking the quark interaction into account. Dashed lines correspond to 
the effective potential V in (5). Bold lines denote the negative-energy part of the quark propagator. 



the subleading order. The vertex functions look like 

r«(p, q) = u f (p f )^(l - 7 5 H(g fe )(27r) 3 <5( P(/ - q 9 ), 

and 



(12) 



e b (k) + e b (p q ) 



7i°V(p 9 - q 9 ) 



+V( Pg - q 



A ( f ~\k') 



qJ e f (k') + e f (q f ) 



TiW 1 -ll)}ub{qb)u q {q q ), (13) 



where the superscripts "(1)" and "(2)" correspond to Figs. 1 and 2, k = pf — A; k' 
q b + A; A = Pf — p Bc ; 



A<">(p) 



e(p) - (m 7 ° + 7 (7P)) 
2e(p) 



Here [20] 



F i=l 

96,, = e b ,g(g)^±]>> W (PB)9\ 

M B i=l 
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and ra W are three four-vectors given by 

" (,) "W = {f-^ + m|tm)}' b = V?^- 

It is important to note that the wave functions entering the weak current matrix element 

(11) are not in the rest frame in general. For example, in the B meson rest frame (ps = 0), 
the final meson is moving with the recoil momentum A. The wave function of the moving 
meson f Fa is connected with the wave function in the rest frame ^fo = by the 
transformation [20] 

*fa( P ) = Df(RYjDy\RYj* F0 (p), (14) 

where R w is the Wigner rotation, L A is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D 1 ^ 2 (R) in spinor representation is given by 

(J ^D^(RYj = S-\p q , c )S(A)S(p), (15) 

where 

is the usual Lorentz transformation matrix of the four-spinor. 

The general structure of the current matrix element (11) is rather complicated, because it 
is necessary to integrate both with respect to d 3 p and d 3 q. The 5-function in the expression 

(12) for the vertex function permits to perform one of these integrations. As a result 
the contribution of to the current matrix element has the usual structure of an overlap 
integral of meson wave functions and can be calculated exactly (without employing any 
expansion) in the whole kinematical range, if the wave functions of the initial and final 
mesons are known. The situation with the contribution r^ 2 ^ is different. Here, instead of a 
5-function, we have a complicated structure, containing the potential of the gg-interaction 
in the meson. It contains also the quark energies e q (p) = \Jm 2 + p 2 , which explicitly depend 
on the the relative quark momentum p. The presence of such dependence does not permit 
one, in the general case, to get rid of one of the integrations in the contribution of r^ 2 ^ to 
the matrix element (11). Therefore, it is necessary to use some additional considerations in 
order to simplify calculations. The main idea is to expand the vertex function T^ 2 \ given 
by (13), in such a way that it will be possible to use the quasipotential equation (1) in 
order to perform one of the integrations in the current matrix element (11) and thus express 
this contribution to the decay matrix element through the usual overlap integral of meson 
wave functions. The realization of this strategy differs for the cases of heavy-to-heavy and 
heavy-to-light transitions. 



IV. SEMILEPTONIC B MESON DECAYS TO D MESONS 

For the description of semileptonic B decays to ground state D mesons (heavy-to-heavy 
transitions) it is convenient to use the HQET parameterization for the decay matrix elements 
[4]: 

( DK)icr^(,)) 

VM D M B +l ; 
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{D(v')\crh 5 \B(v)) = 0, 



(D*(v',e)\c^b\B(v)) 
\J Mp* Mb 
(D*(v',e)\cYl 5 b\B(v)) 
\J Mo* Mb 



h Al (w + - (h A y + h A3 v'»)(e* ■ v), 



(16) 



where v (v') is the four- velocity of B (D^) meson, e M is the polarization vector of the 
final vector meson, and the form factors hi are dimensionless functions of the product of 
four-velocities 



w — v • v 



, M l + Km - v 2 



2M B M D («) 

and q = Pb — PdM is the momentum transfer from the parent to daughter meson. 
In HQET these form factors up to l/rnq are given by [4] 

h + =£ + (e c + e b ) [2 X i - 4(iu - l)x 2 + 12 X3 ] , 

h_=(e c -e b ) [2^3 -M] , 



h v =i + e c 2 X i - 4x3 + M 



h Al =£ + e c 



2xi - 4xs + 



w + 1 



1 _ 



2xi - A(w - 1)X2 + 12X3 + M ~ 2^3 



+ e b 



2xi - A(w - 1)X2 + 12X3 + ^— J (M ~ 26) 

w + 1 v 7 



4X2 



w + 



T (A« + &) 



+e b 



w 



(n) 



2X1-4X2-4X3 + 

w+1 w+1 
2xi - 4(w - l)x 2 + 12X3 + M - 26 
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(17) 



The calculation of the weak transition B ^ matrix elements using the heavy quark 
expansion shows that all model independent HQET relations are satisfied in our model. In 
the limit of an infinitely heavy quark all form factors are expressed through the Isgur-Wise 
function [2] 



h+(w) = h Al (w) = h Aa (w) = h v (w) = £(w) 
h-(w) = h A2 (w) = 0. 

In our model this function is given as the overlap integral of meson wave functions [8] 



w + 



IW 



w + 1 



e A *b(p), 



(18) 



(19) 



where = A/ V A 2 is the unit vector in the direction of A = Mp v' — Mgv. In order to 
fulfill the HQET relations (17) in the first order of the heavy quark I/tuq expansion it is 
necessary to set (1 — e)(l + k) = 0, which leads to the vanishing long-range chromomagnetic 
interaction. This condition is satisfied by our choice of the anomalous chromomagnetic quark 
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FIG. 3: The Isgur-Wise function £(w). 

moment k — — 1. To fulfill the HQET relations at second order in I/tuq it is necessary to 
set £ = — 1 [8]. This gives an additional justification, based on the heavy quark symmetry 
and heavy quark expansion in QCD, of the choice of the characteristic parameters of our 
model. In the infinitely heavy quark mass limit the wave functions of initial and final 
\&£) heavy mesons coincide. As the result the HQET normalization condition [4] 

= i 

is exactly reproduced. The first order Isgur-Wise functions are given by [8] 



X2(w) 

XsO) 



(A - m q ) 1 + - 



3w + 1 



A— ««), 



w + 1 
1 A 



32 w 

1 T w ■ 
— A- 



en, 



(20) 



where the HQET parameter A 
heavy meson 



16 w + 1 

M — rriQ is equal to the mean energy of a light quark in a 
A = (e q ) ~ 0.56 GeV. 

The functions Xi and X3 explicitly satisfy normalization conditions at the zero recoil point 
[21] 

Xl (l) = X3 (l)=0, 

arising from vector current conservation. 

In Figs. 3-5 we plot the Isgur-Wise functions calculated with numerical wave functions 
determined in the process of their mass calculations [11]. Near the zero recoil point of the 
final meson w = 1 the Isgur-Wise function can be written as 



p 2 (w - 1) +c(w - l) 2 + 



(21) 
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FIG. 4: The first order functions £s(w)/A (bold line) and xi(w)/A (solid line). 
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FIG. 5: The first order functions — X2(w)/A (bold line) and X3( w )/^ (solid line). 



where p 2 = — [d£(w) /dw] w= i ~ 1.04 is the slope and c = (1 / '2)[d 2 ^(w) / 'd 2 w] w= i ~ 1.36 is the 
curvature of the Isgur-Wise function. The slope p 2 can be compared to the recent quenched 
lattice QCD evaluation [22]: p 2 = 0. 83i 11 i 2 2- Note that both the slope and curvature of the 
calculated Isgur-Wise function satisfy all known lower bounds (see Ref. [23] and references 
therein) . 

The differential semileptonic decay rate B — > Dlv for the massless leptons is given by 

[4] 



dT 



G 2 



\V cb \ 2 Ml (w 2 - lf'\M B + M D f F 2 ( 



w 



dw 487T 3 

where Gp is the Fermi constant, and the form factor Fd(w) is defined by 



Fn(w) 



h+(w) — 



1 — r 
1 + r 



M. 



D 



M, 



B 



(22) 



(23) 
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FIG. 6: Comparison of experimental data and predictions of our model for the product Fi)(w)\V c b\. 
Dots represent CLEO data [24] and diamonds show Belle data [25]. Solid lines are predictions of 
our model for \V cb \ = 0.044, 0.039, 0.034 (from top to bottom, respectively). 



Near the zero recoil point the form factor Fd(w) has the following expansion 

F D {w) = F D {l)[l - p D (w - 1) + c D (w - l) 2 + • • ■], (24) 

where the value of F D (1), calculated by using the unexpanded in l/mg expressions (Al)- 
(A3) of Ref . [9] , is equal to 

F D (1) = 0.966. (25) 

Quenched lattice QCD calculations [26] give the value F D (1) = 1.058 ± 0.016 ± 0.03^ ^, 
while Ref. [27] predicts F D (1) = 0.98 ± 0.07. 

The slope of the form factor Fd(w) at zero recoil w — 1 in our model 



2 1 dF D (w) 

Pd ~ 



Fd(w) dw 



0.88, (26) 



w=l 



is in good agreement with experimental values p 2 D = 0.76 ± 0.16 ± 0.08 [24] and 
p 2 D = 0.69 ± 0.14 [25], obtained by using the linear fit of the data. The curvature of the 
form factor F D (w) is equal to c D = [1/ 1 (2F D (l))][d 2 F D (w) / 'd 2 w] w=1 ~ 0.75. 

In Fig. 6 we compare the results of our model for the product Fo{w)\V c b\ with the recent 
experimental data of CLEO [24] and Belle [25]. It is seen that the form factor -Fd(w) 
dependence on w in our model is in good agreement with measurements. The combined fit 
of CLEO (F D (l)\V cb \ = 0.039 ± 0.002) and Belle (F D (l)\V cb \ = 0.041 ± 0.003) data leads to 
the value of the product of the form factor and CKM matrix element 

F D (l)\V cb \ = 0.040 ±0.002. 
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Using our prediction (25) for the form factor F D (1) we find the value of CKM matrix element 



\V cb \ = 0.0415 ±0.0020. 
The differential semileptonic decay rate B — > D* I v is defined by [4] 



dw 



48tt 3 



\V cb \ 2 (M B - M D .y M 6 D , J(w 2 -l)(w + 1) 



x 



1 + 



Aw 1 — 2wr* + r* 2 



w + 1 (1 — r 
where the form factor Fp* (w) is given by 



*\2 



F 2 », 



M 



D* 



F D » (w) = h Al (w) 



H\{w) + H 2 _(w) + H$(w) 



1 + 



Aw 1 — 2wr* + r 



*2 • 



^ w + 1 (l-r*y 



The helicity amplitudes Hj(w) 



H±(w) = 



VT — 2wr* + r* 2 



tw + 1 



tfoH = l + ^4[l-i? 2 H]r 



are expressed through form factor ratios 



h v (w) 
h Al (w) : 

h A3 (w)+r* h A2 



h Al (w) 



(27) 



(28) 



(29) 



(30) 



(31) 



In the limit mq — > oo, i?i = R 2 — 1 due to spin-flavour symmetry [4]. Taking into account 
of l/rriQ corrections breaks down this symmetry relation. In our model using unexpanded in 
1 / rriQ formulas for the form factors [9] we get the following expressions near the zero recoil 
point w — 1 



0.86(w- l) + 0.72(w- 1) 2 + - 
-0.66(w-l) + 0.17(w-l) 2 + 



h Al (w) = 0.918[1 • 
F D .(w) = 0.918[1 
i?i(w) = 1.39 - 0.23(w - 1) + 0.21(w - l) 2 + 
R 2 ( w ) = 0.92 + 0.12(iu - 1) - 0.07(w - l) 2 + 



(32) 



It is necessary to note that the behaviour of the form factor h Al (w) predicted by our model 
agrees, in general, with the parameterization [28] in the whole kinematical range. On the 
other hand, the expansion of the form factor ratios R\ and R 2 is close to the QCD sum rule 
results [4]. Lattice QCD calculation [29] gives the value of h Al (l) = 0.9130j#oira3l:o302- Out 
result for the difference of the slope parameters 



PM ~ Pi* = 0.20 
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w 



FIG. 7: Comparison of experimental data and predictions of our model for the product Fjj* (w) \ V c b \ ■ 
Dots show CLEO data for B+ -» D*°l~u, triangles - CLEO data for B° -» D*+l~v [31], diamonds 
- Belle data [32], squares - BaBar data [33]. Solid lines show predictions of our model for \V c b\ = 
0.044, 0.039, 0.034 (from top to bottom, respectively). 



is in agreement with previous calculations 0.21 [27] and 0.17 [30]. Comparing Eqs. (26) and 
(32) we find the value of the slope difference 

Pm -pI = -°- 02 > 

which coincides with the one found in Ref. [30]. 

In Table I we give our predictions for i?i(l), -^2(1) and for the slope of the form factor 
h Al and available experimental data [31-34]. Our results for R\ and R 2 agree well with data. 
The calculated slope p\ is within experimental error bars for the values obtained by using a 
linear fit. Comparison of our predictions for the product Fjj*(w)\V c i,\ with the experimental 
data from CLEO [31], Belle [32] and BaBar [33] is given in Fig. 7. In general there is good 
agreement between the calculated form factor behaviour and available experimental data. 
The values of the product Fn*{l)\V c b\, obtained by using our form factors are compared in 
Table I with the ones based on the parameterization of the form factor Ha x from Ref. [28]. 
Our model leads to the values of the product Fd*(1)|14&| approximately 10-15% lower, than 
the values obtained using the parameterization [28] . A combined fit of all above mentioned 
experimental data in the framework of our model gives 

F D *(l)\V cb \ = 0.0343 ±0.012. 

Using our value of F D *(1) = h Al (l) = 0.918 (see (32)) we get 



(33) 
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TABLE I: Comparison of our model predictions for the ratios -Ri(l), B,2(l), the slope of the form 
factor fiAi for w = 1 and values of the product Fr>*(l)\V cb \ with experimental data. 



our 



CLEO [31] 



BaBar [33] 



Belle [32] 



DELPHI [34] 



Ri 
R2 



Pi 



1.39 
0.92 

0.86 



F D *\V cb \ 0.0343(12) 



1.18(30)(12) 
0.71(22)(7) 

f 0.91(15)(6) a 

\ 1.61(9)(21) 6 

0.0360(20) c 

0.0431(13)(18) 6 



1.396(60)(44) 

0.885(40)(26) 

0.79(6) a 

1.145(59)(46) 6 

0.0328(5) c 

0.0376(3)(16) 6 



J 0.81(12) a 
\ 1.35(17) 6 

0.0315(12) c 

0.0354(19)(18) 6 



1.39(10)(33) t 



0.0377(11)(19) 6 



a linear fit of experimental data. 

b fit using the form factor Ha 1 parameterization [28] . 
c fit using form factor predictions of our model. 



TABLE II: Comparison of theoretical predictions for the ratios R\(l), #2(1) and their derivatives 
^i(l), ^(1)- 



Ref. R[{1) R 2 (l) i%(l) 

~om L39 -0.23 092 0.12 

[30] 1.25 -0.10 0.81 0.11 

[27] 1.27 -0.12 0.80 0.11 

[4] 1.35 -0.22 0.79 0.15 

[35] 1.15 0.94 

[36] 1.01(2) 1.04(1) 



The theoretical results for the ratios -Ri(l), #2(1) an d their derivatives R[(l), R^iX) are 
confronted in Table II. In general there is a reasonable agreement both for the form factor 
ratios and their slopes. Our values are very close to the ones found using HQET and QCD 
sum rules [4]. 

By integrating the expressions for the differential decay rates (22), (28), we get predictions 
for the total decay rates in our model 

T(B -> Dlu) = 9A8\V ch \ 2 ps" 1 , 

T(B -> D*lv) = 24.9\V cb \ 2 ps" 1 . (34) 

Taking mean values of lifetimes [37]: t b o = 1.530 x 1(T 12 s and r B + = 1.671 x 10~ 12 s, we 
find 



BR(B° - 


■+D + r 




= 14.5] 


v cb \ 2 


BR(B + - 






= 15.8; 


v cb \ 2 


BR(B° - 


► D* + r 




= 38.0; 


v cb \ 2 


BR(B + - 






= 41.3: 


v cb \ 2 



The comparison of theoretical and experimental branching ratios [37] leads to the following 
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values of the CKM matrix element \V cb \: 

BR(B° -> D + rv) cxp = 0.0214 ± 0.0020 \V cb \ = 0.038 ± 0.002, 

BR(B + -> D°l + u) cxp = 0.0215 ± 0.0022 jv^j = 0.037 ± 0.002, 

BR(B° -> £>* + Z - i/) exp = 0.0544 ± 0.0023 jv^j = 0.038 ± 0.001, 

£/?(£+ -> D*°l + u) cxp = 0.065 ± 0.005 jv^j = 0.040 ± 0.002, (36) 

which are in good agreement with each other and with values (27), (33), found from the 
form factor analysis. Thus the averaged \ V c b\ over all presented experimental measurements 
of semileptonic decays B — > Dev and B — > £>*ez/ is equal to 

\V ch \ = 0.0385 ±0.0015 (37) 

in good agreement with PDG [37] 

\V cb \ = 0.0409 ± 0.0018 (exclusive). 



V. SEMILEPTONIC B DECAYS TO LIGHT MESONS 



The matrix elements of weak current J w governing the weak B decays to the light 
pseudoscalar meson (P = n) is parameterized by two invariant form factors. It is convenient 
to use the following decomposition: 



(P(p F )\q^b\B(p B )) = f + (q 2 ) 



Pb+Pf 



Ml 



+ H^ 



Ml „ 

— ^, (38) 



where q = Pb — Pf] Mb is the B meson and Mp is the final pseudoscalar meson mass. 

The corresponding matrix elements for the weak B decays to the light vector meson 
(V = p) can be parameterized by four form factors: 



{V(p F )\qYb\B(p B )) 



2iV(q 2 
Mb + Mi 



v 



-e^elpBpPFa, 



(39) 



(V(p F )\qr^b\B(p B )) = 2M v A (q 2 ) e —^ 

e* • q 



q» + {M B + M v )A l {q 2 )[e^ 



-A 2 (q 2 ) 



Mb + My 



Pb+Pf ~ 



M 2 B 



(40) 



where My and are the mass and polarization vector of the final vector meson. At the 
maximum recoil point (q 2 = 0) these form factors satisfy the relations: 



A o (0) 



/+(0) = /o(0), 

Mr + Mi/ Mr — Mi 

-^i(0)- \, — A 2 (0). 



■2 My V ' 2 My 

For massless leptons form factors fo and A$ do not contribute to the semileptonic decay 
rates. However they give contributions to the nonleptonic decay rates in the factorization 
approximation. Note that this parameterization is completely equivalent to the one of HQET 
(16), and corresponding form factors can be easily expressed through each other. 
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In this section we calculate semileptonic decay rates of the heavy B meson into light 
meson, B — > 7i(p)eu. The final meson in these decays contains light quarks (u, d, s) only, 
thus in contrast to decays to heavy D mesons, considered previously, the application of the 
expansion in inverse powers of the final active quark is not justified. The calculation of the 
contribution of the vertex function (12) to the decay matrix element of the weak current 
(11) can, as it was already in detail discussed above, be carried out exactly, due to the pres- 
ence of 5-function, and does not require any expansion. The calculation of the contribution 
I^ 2 ) is significantly more difficult, since the expansion only in inverse powers of the heavy 
6-quark mass from the initial B meson retains the dependence on the relative momentum in 
the energy of the final light quark. Such dependence does not allow one to perform one of 
the integrals in the decay matrix element (11) using the quasipotential equation. However 
the final light meson has a large (compared to its mass) recoil momentum (A = — p#, 
I A max | = (M§ — Mp)/{2M B ) = M B /2 ~ 2.6 ) almost in the whole kinematical range except 
the small region near q 2 = q^ x (|A| = 0). This also means that the recoil momentum of 
the final meson is large with respect to the mean relative quark momentum |p| in the meson 
(~ 0.5 GeV). Thus one can neglect |p| compared to |A| in the final light quark energy 
e q (p + A) = \Jm 2 + (p + A) 2 , replacing it by e g (A) = yVn 2 + A 2 in expressions for the T^ 2 ) 
contribution in accord with the large-energy expansion (see Introduction). This replace- 
ment removes the relative momentum dependence in the energy of the light quark and thus 
permits to perform one of the integrations in the T^ 2 ^ contribution using the quasipotential 
equation. Note that the relatively small value of this contribution, related to its proportion- 
ality to the binding energy in the meson, and its predictable momentum dependence allow us 
to extrapolate it to the whole kinematical range. The numerical analysis (see below) shows 
that such extrapolation induces insignificant uncertainties in the final results for decay rates. 

It is convenient to consider semileptonic decays B — > {it,p)ev in the B meson rest frame. 
Then calculating decay matrix elements it is necessary to take into account the relativistic 
transformation (14) of the final meson wave function from the rest frame to the moving 
one with the momentum A. Applying the method described above, we find expressions for 
the decay matrix element and determine the corresponding form factors. They have the 
following structure: 

(a) B — > 7r transitions 

U(q 2 ) = /iV) + e/fV) + (1 - e)f^ (g 2 ), (41) 

fo(q 2 ) = rfV) + £/o 5(2 V) + (1 " s)f^\q% (42) 

(b) B — > p transitions 

V(q 2 ) = V^(q 2 )+eV s ^(q 2 ) + (1 - e) V v ^ (g 2 ), (43) 

Mq 2 ) = A^\q 2 ) + eA s ^\q 2 ) + (1 -e)< (2 V), (44) 

Mq 2 ) = A?(q 2 )+eA^\q 2 ) + (1 - e)A^ 2 \q 2 ), (45) 

Mq 2 ) = A^(q 2 ) + eA^ 2 \q 2 ) + (1 - e)< (2 V), (46) 

where expressions for f^ , f+p^, A^ 12 , Aq'J^, and V s,v ^ are rather cumbersome 
and can be found in the Appendix to Ref. [10]. The subscripts "(1)" and "(2)" correspond 
to Figs. 1 and 2, S and V denote scalar and vector potentials of the gg-interaction. Let us 
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remind that the mixing coefficient e of vector and scalar confining potentials is equal to —1 
in our model. Note that form factors (41)-(46) in the limit of the infinitely heavy b quark 
mass and large recoil of the final light meson explicitly satisfy all symmetry relations [6, 7] 
imposed by the large energy effective theory. 

In order to increase the precision and reliability of our calculations compared to our 
previous consideration [10] we do not perform a further expansion of form factors in inverse 
powers of the heavy quark mass. Moreover, the q 2 dependence of form factors is explicitly 
determined by these formulas and thus no ad hoc parameterization is necessary. We also 
use the numerical wave functions found in the meson mass spectrum calculations [11, 12] 
instead of trial (Gaussian) wave functions used in Ref. [10]. To check the precision of the 
extrapolation of the form factors in the region of small recoil of the final light meson (A = 0, 
q 2 = g^ax) we perform the additional consideration of form factors in this region. In this 



analysis the simplifying substitution of the light quark energy e q (p) = ^Jp 2 + m 2 q by the 
on-shell center of mass energy E q = (M 2 — m 2 Q + m 2 )/(2M) (3) in the contribution of the 
vertex function is used. Such replacement is valid in the small region around the zero 
recoil point. As a result in this point one gets expressions for the form factors similar to 
the formulas given in the Appendix of Ref. [38] with obvious substitutions. Such calculation 
showed that the obtained values of form factors for semileptonic decays B — > (ir, p)ev as 
well as their slopes at zero recoil q 2 = g^ ax are in good agreement with the results found 
from the extrapolation of Eqs. (41)-(46). The deviations are less than 1% confirming the 
reliability of such extrapolation and of final results which are based on it. 

The semileptonic B — > [n,p)ev decay from factors in our model can be approximated 
with good accuracy by the following expressions [35]: 

^)F(q 2 ) = U(q 2 ),V(q 2 ),A (q 2 ) 

*V) = ~ „ 2N 7 (47) 




1 - 0-1T79 - + a 2- 



M 2 B , M%, 



(b) F(q 2 )=A 1 (q 2 ),A 2 (q 2 ) 



/'V) = - -4^ *TT> ( 48 ) 



q- q- 



where M = Mb for Aq and M = Mb* for all other form factors; the values F(0) and <j\^ 
are given in Table III. The difference of fitted form factors from the calculated ones does 
not exceed 1%. 

In Table IV we give a comparison of the predictions for the form factors of semileptonic 
decays B — > (jr, p)ev at maximum recoil point q 2 = 0, calculated in our model with the 
results of light cone QCD sum rules (LCSR) [39-41], the quark model (QM) [35], using 
relativistic dispersion relations and two recent lattice QCD (LQCD) calculations [42, 43] with 
dynamical light quarks. Note that our new results for form factors at this point coincide 
with previous calculations [10] within errors caused by the Gaussian parameterization of 
the heavy-light meson wave functions used in Ref. [10]. In Ref. [44] a model independent 
constraint for the product |V^&|/+(0) = (7.2 ± 1.8) x 10~ 4 was obtained using the soft- 
collinear effective theory and B — > nn data, which for their final value of \ V u b\ yields /+(0) = 
0.227 ±0.047. 
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TABLE III: From factors of semileptonic decays B — > (ir, p)ev calculated in our model. Form 
factors f+(q 2 ), V(q 2 ), A (q 2 ) are fitted by Eq. (47), and form factors A 1 (q 2 ), A 2 (q 2 ) are fitted by 
Eq. (48). 



B -► 7T S -► p 





/+ 


fo 


V 


A 


A! 


A 2 


F(0) 


0.217 


0.217 


0.295 


0.231 


0.269 


0.282 




0.378 


-0.501 


0.875 


0.796 


0.54 


1.34 


0"2 


-0.41 


-1.50 





-0.055 





0.21 



TABLE IV: Comparison of theoretical predictions for the form factors of semileptonic decays 
B — > (7r, p)ezv at maximum recoil point g 2 = 0. 





/+(0) 


V(0) 


4>(0) 


A 1 {0) 


A 2 (0) 


our 


0.217 


0.295 


0.231 


0.269 


0.282 


LCSR [39, 40] 


0.258(31) 


0.323(29) 


0.303(28) 


0.242(24) 


0.221(23) 


LCSR [41] 


0.25(5) 


0.32(10) 




0.24(8) 


0.21(9) 


QM [35] 


0.29 


0.31 


0.29 


0.26 


0.24 


QM [10] 


0.20(2) 


0.29(3) 




0.26(3) 


0.31(3) 


LQCD(FNAL)[42] 


0.23(2) 










LQCD(HPQCD)[43] 


0.27(5) 











The q 2 form factor dependence is plotted in Figs. 8 and 9. In Fig. 8 we also show recent 
lattice results for the form factors f+(q 2 ) and fo{q 2 ) [42, 43]. It is clearly seen from this 
figure that the behaviour of the form factor f+(q 2 ) agrees with lattice computations within 
errors, while our form factor fo(q 2 ) lies somewhat lower than lattice data. In this figure we 
also show the LCSR result for the value of form factors at q 2 = [39]. It agrees with our 
model prediction within errors. 

We can also check the consistency of the obtained q 2 behaviour of the form factor / + by 
comparing its calculated value at q^^ with model independent results of chiral perturbation 
theory (ChPT). For the pion recoil energy E n ~ m n ChPT predicts [44] 



gf B M E 



2f n (E„ + Mb* - M B ) 



(49) 



where g ~ 0.5 is B*Bir coupling [44], the decay constant fs is equal 189 MeV in our model 
[45]. The first corrections scale as E w /A, where A ~ 600 MeV is the mass splitting to the 
first radially excited l~ state above the B*. Substituting these values one gets the following 
prediction [44] 

/+(&*) = 10.38 ±3.63, (50) 

which is in agreement with our model result f + (q 2 iax ) = 10.9. 

On the other hand, the form factor /o in the soft-pion limit p — > and m 2 — > is related 
to the ratio of the B and ir decay constants [39, 46] 



f ( 2 \ f B 

Jn 



(51) 
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FIG. 8: Form factors of the semileptonic decay B — * ixev in comparison with lattice calculations: 
f+(q 2 ) is given by the upper plot (our), squares (FNAL)[42], diamonds (HPQCD)[43]; fo{q 2 ) is 
given by the lower plot (our), stars (FNAL)[42], triangles (HPQCD)[43]. The LCSR value for from 
factors at q 2 = is plotted by a circle [39]. 



7 2 ) 

imax/ 



1.45 



This relation with the above values of the decay constants gives the result fo(q 2 
again in good agreement with the prediction of our model /o(O'max) = 1-36. 

The differential semileptonic decay rates can be expressed in terms of these form factors 

by 

(a) B — ► Pev decay (P = n) 



r<2 A3IT/ 12 



(b) B -> Veu decay (V = p) 



AV AIT/ I 2 a 2 

where Gp is the Fermi constant, V q b is CKM matrix element (q = u), 

A = |A 



\ 



(M 2 + M 2 - g 2 ) 2 



4Mi 



M 2 y . 



Helicity amplitudes are given by the following expressions 



H ± (q 2 



2M B A 
Mb + Mv 



(M B + M v 
2M B A 



H (q 2 



(M B + M y )(Ml - My - g 2 )A(g 2 



4M 2 A 2 , , 2N 



Mb + My 



(52) 



(53) 



(54) 
(55) 
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FIG. 9: Form factors of the semileptonic decay B — ► pew. V(q 2 ) is plotted by the solid line, A\(q 2 
- by the bold line, A2(q 2 ) _ by dashed line, and A${q 2 ) - by long-dashed line. 



The decays rate in transversally and longitudinally polarized vector mesons are defined 



by 



dr 7 



dr, 



+ 



dr L 

dg 2 

dr_ 



G 2 F A\V n 



, 2 a 2 
qb\ Q 



96tt 3 



M 2 



B 



G 2 F A\V qb \ 2 q 2 



dq 2 dq 



967T" 



M 2 B 



\H*{q 2 )? 
(\H + {q 



2\|2 



+ \H-{q 2 )? 



(56) 
(57) 



dq 2 

Integration over q 2 of these formulas gives the total rate of the corresponding semileptonic 
decay. These rates are presented in Table V. In this Table we also give the values of 
the partial decay rates integrated over two intervals q 2 < 16 GeV 2 and q 2 > 16 GeV 2 in 
comparison with the evaluations of LCSR in the first interval and LQCD in the second 
interval [47] . This is related to the fact that the predictions of these approaches are reliable 
only in the above mentioned regions. The results presented in Table V show that our 
calculations agree well with the lattice evaluations, while our values are somewhat lower 
than LCSR ones which have relatively large errors. Note that our model gives the ratio of 
semileptonic B decay rates into p meson with longitudinal polarization to the corresponding 
rate with transverse polarization Yl/Yt = 0.46. 

Using mean lifetimes of B mesons [37]: tb° = 1.530 x 10~ 12 s and = 1.671 x 10~ 12 s, 
we find 



BR(B° - 




= 8.34|K&| 2 


BR(B + - 




= 4.47|K fe | 2 


BR(B° - 


- p+l-u) 


= 20.1\V ub \ 2 


BR(B + - 


- p°l + u) 


= l0.7\V ub \ 2 



(58) 

Comparison of these predictions with experimental data [37] leads to the following values of 
\V ub \: 



BR(B° -> 7r + rz/) cxp = (1.36 ± 0.15) x 10" 



| Kb | = (4.04 ±0.25) x 10" 
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TABLE V: Comparison of theoretical predictions for the rates of semileptonic decays B — > (tt, p)eu 
(in iKfePps- 1 )- 









r 




]V<16 GeV 2 




pg 2 >16 GeV 2 


Decay 


our 


LCSR 


FNAL 


HPQCD 


our LCSR 


our 


FNAL HPQCD 






[39] 


[42] 


[43] 


[39] 




[42] [43] 


B — > 7rez^ 


5.45 


7.74(2.32) 


6.24(2.12) 


6.03(1.94) 


3.68 5.44(1.43) 


1.77 


1.83(50) 1.46(35) 


£> — > pez^ 


13.1 








10.5 


2.60 





BR{B + -> 7r°/ + z/) cxp = (7.4 ± 1.1) x 10" 4 \V uh \ = (4.07 ± 0.30) x 1(T 3 , 

BR{B° -> p + rv) cxp = (2.3 ± 0.4) x 10" 4 \V ub \ = (3.38 ± 0.30) x 10" 3 , 

BR{B + -> p°/+z/) cxp = (1.24 ± 0.23) x 10" 4 \V ub \ = (3.39 ± 0.33) x 10" 3 . (59) 

Decays of the neutral and charged B mesons give very close results for \V u i\, while the 
averaged values of \V u f,\, extracted form the decay B — > pev are approximately 16% lower 
than corresponding values, found from decay B — > -nev. Note that the recent CLEO [48] 
measurement of the decay branching ratio BR(B° — > p + l~v) cxp = (2.91 ± 0.54) x 10~ 4 gives 
\V u b\ = (3.81 ± 0.35) x 10~ 3 which is close to the one extracted from the B — > nev decay. 

Recently significant progress has been achieved in the experimental determination of the 
differential decay B — > (jr, p)ev rate dependence on q 2 . CLEO [48, 49], BaBar [50-52], Belle 
[53] measured partial decay rates in relatively narrow intervals of q 2 . In Figs. 10 and 11 we 
present the comparison of our model predictions for partial branching ratios of B — > (tt, p)ev 
decays with experimental data. Plotting histograms in Figs. 10 (a)- (c) and Figs. ll(a),(b) 
we used the value of \V u b\, extracted from the total rate in the corresponding experiments. 
It is clearly seen that our predictions agree well with almost all data for decays of neutral as 
well as charged B mesons. Using these experimental data on partial and total semileptonic 
B — > {ji,p)eu decay rates it is possible to extract averaged values of \V u b\: 

B^ixev | Kb | = (4.02 ±0.10) x 10~ 3 , 

B -> pev IK&I = (3.33 ± 0.20) x 10" 3 , (60) 

which are in good agreement with the ones, extracted from averaged total decay rates (59). 

Finally, averaging over all above mentioned experimental data we get the value of |T4b| 
in our model 

|Kb| = (3.82 ±0.20) x 10~ 3 (61) 
in good agreement with PDG [37] 

|K&| = (3.84+°;^) x 10~ 3 (exclusive). 

VI. CONCLUSION 

We calculated weak decay form factors and decay rates of different semileptonic B decays. 
Decays both into heavy and light 7r(p) mesons were considered. 

First, it was shown that our relativistic quark model gives a reliable description of the 
heavy-to-heavy semileptonic transitions B — > D^lv. All model independent HQET rela- 
tions are reproduced. The model allows to express corresponding leading and subleading 
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FIG. 10: Comparison of theoretical predictions for partial decay rates with experimental data: 
(a)-(c) (l/F tot )(J Aq2 {dr(B -> Trev) /dq 2 )dq 2 ) x 10 4 from Refs. [48, 50, 53]; (d), (e) ABR{B -> 
nev) / BR{B — > nev) from Refs. [51, 52], respectively. Lower histograms on (a), (b) and diamonds 
are theoretical and experimental values for decays of charged B meson B + — » 7r°e + zA All other 
histograms and squares are theoretical predictions and experimental data for decays of neutral B 
meson B° — ► 7r + e~zA 



order Isgur-Wise functions through overlap integrals of meson wave functions. These wave 
functions were determined previously in the process of the meson mass spectrum calculations. 
From the comparison with the experiment and predictions of other theoretical approaches 
it follows that our model correctly reproduces the q 2 behaviour of form factors. Calculated 
decay rates and branching fractions are in good agreement with data and give very close 



23 



1 . 6 
1 . 2 

. 8 

. 4 




(a) 



10 



1 . 6 
1 . 2 

. 8 

. 4 



15 



20 




(b) 



10 



15 



20 




FIG. 11: Comparison of theoretical predictions for partial decay rates with experimental data: (a), 
(b) {l/F tot )(J Aq2 {dr(B -► pev)/dq 2 )dq 2 ) x 10 4 from Refs. [49, 53]; (c) ABR(B -► pev)/BR{B -> 
pe^) from Ref. [51], respectively. Lower histograms on (b) and diamonds are theoretical and 
experimental values for decays of charged B meson B + — > p e + u. All other histograms and squares 
are theoretical predictions and experimental data for decays of neutral B meson B° — ► p + e~v. 



values of the CKM matrix element |V^,| extracted from different decay processes. 

Secondly, the form factors of the heavy-to-light semileptonic B decays were calculated. 
The consideration was done with the systematic treatment of all relativistic effects, which 
are very important for such transitions. Particular attention was paid to the inclusion 
of negative-energy contributions and to the relativistic transformation of the meson wave 
function from the rest to the moving reference frame. The q 2 dependence of the form 
factors was explicitly determined without using any ad hoc assumptions. The validity of the 
form factor extrapolation, which is necessary only within the small region near the point 
of zero recoil of the final light meson, was crosschecked by an additional calculation of the 
form factor values in this particular point. The decay form factors are again given by the 
overlap integrals of the B and tt, p meson wave functions, which are know from the previous 
calculations of the meson mass spectra. The q 2 behaviour of the form factors is in agreement 
with both unquenched lattice QCD calculations and predictions of light cone QCD sum rules 
within the ranges of the validity of these approaches. All this significantly improves the 
reliability of the obtained results. The comprehensive comparison of the predictions with 
recent experimental data both on total and partial decay rates allowed the extractions of 
the CKM matrix element \V u b\ which values are rather close in different decay channels. 
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The evaluation of the theoretical uncertainties represents an important problem. Unfor- 
tunately, it is not easy to estimate them in the framework of the adopted approach. The 
theoretical errors within the model, which come from the neglected higher order terms in 
the 1/m.Q expansion for heavy-to- heavy transitions and from the form factor extrapolation 
in the region of zero recoil for heavy-to-light transitions, can be easily estimated and are 
less than 3%. The main difficulty is related to the uncertainty of the quark model itself. 
However, our previous experience in calculating a vast set of different properties of hadrons 
within our model indicates that such uncertainties should not exceed 5%. Therefore, adding 
these errors in quadrature we find for the CKM matrix elements the following final values 
in our model: 

\V cb \ = (3.85 ± 0.15 ± 0.20) x 10" 2 , 

\V ub \ = (3.82 ± 0.20 ± 0.20) x 10^ 3 , (62) 
where the first error is experimental and the second one is theoretical. 
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